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Fluctuation of Coulomb blockade peak spacings in large two-dimensional 
semiconductor quantum dots are studied within a model based on the elec- 
trostatics of several electron islands among which there are random inductive 
and capacitive couplings. Each island can accommodate electrons on quan- 
tum orbitals whose energies depend also on an external magnetic field. In 
contrast with a single island quantum dot, where the spacing distribution is 
close to Gaussian, here the distribution has a peak at small spacing value. 
The fluctuations are mainly due to charging effects. The model can explain 
the occasional occurrence of couples or even triples of closely spaced Coulomb 
blockade peaks, as well as the qualitative behavior of peak positions with the 
applied magnetic field. 
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I. INTRODUCTION 



Recently, it became apparent that the physics exposed in the addition spectra of quan- 
tum dots is rather rich, and hence its investigation is at the focus of both experimental and 
theoretical studies. The present work concentrates on the distribution of spacings between 
Coulomb blockade peaks in large semiconductor quantum dots. Coulomb blockade is evi- 
dently one of the hallmarks of mesoscopic physics. The experimental achievement of tracing 
an addition of a single electron to a quantum dot and the appearance of isolated conductance 
peaks led to the concept of single electron transistors. After the origin of Coulomb blockade 
peaks has been elucidated, investigation is directed toward more subtle questions like their 
heights, widths and spacings. The underlying physics is related to the ground state energy, 
chemical potential and inverse compressibility of a few electron islands coupled capacitively 
to its environment, as well as fluctuations of these quantities with the number N of electrons 
on the dot. 

As far as the distribution of spacings between adjacent Coulomb blockade peaks is concerned, 
the question can be stated as follows: According to the simplest picture in which the quantum 
dot is regarded as a single electron island whose coupling with the leads is through its 
capacitance C, the total potential energy of a quantum dot is Q 2 /2C — V g Q, where V g is the 
corresponding gate voltage. The conductance peaks occur at those values of V g for which 
CV g = e(N + 1/2), where e is the electron charge (henceforth e = —1). For this value of V g 
the addition of an electron to the dot (which contains N electrons) does not cost any charging 
energy. The position of the N th Coulomb blockade peak is then a linear function of N, and 
therefore the spacing should be a constant 1/C, independent of N. Recent experiments |l],[2[ 
indicate however that spacing between Coulomb blockade peaks in small quantum dots is in 
general not constant but, rather, a fluctuating quantity close to Gaussian. The average of 
its distribution approximately coincides with the constant value mentioned above, but the 
elucidation of its standard deviation is still under investigation 

The situation is even more intriguing if the quantum dot is very large. As indicated in a series 
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of recent experiments, the spacing occasionally vanishes, namely, two peaks (and sometimes 
even three peaks) coincide. Moreover, the evolution of peak positions and spacings with 
an applied magnetic field indicates the existence of strong correlations between them 0-0 . 
These observations motivated numerous theoretical models based on the concept of pair 
tunneling || or that of two-electron bound-states in depleted electron islands 0. 
In the present work we examine the scenario according to which a large quantum dot, like 
the one used in the last experiment 0], is, in fact, composed of several electron islands which 
are coupled capacitively among themselves as well as to the leads. Electrons are added in 
such a way that the total potential energy of the dot is minimum. This simple generalization 
of the single island picture leads to a remarkable change in the spacing distribution from 
a Gaussian |],fi|,[| centered around a finite average to one which is large at small spacings. 
When the coupling between islands is zero, the distribution has indeed a maximum at zero 
spacing. If the coupling is non-zero but small, the maximum is slightly shifted away from 
zero, yet leaving an appreciable tail down to zero. These results are short of explaining 
the perfect overlap of peaks, since it requires a delta function component at zero spacing. 
Yet, it leads to the occurrence of couples and sometimes triples of closely spaced peaks, 
similar to the experimental observation. Moreover, the evolution of the peak positions with 
the magnetic field is qualitatively similar to the experimental one. On the other hand, the 
present model does not predict a definite periodicity in the bunching of Coulomb blockade 
peaks with electron number N. In the next section the model is explained, and the results 
of the calculations are presented in section 3. 



II. FORMALISM 

We consider a large isolated two-dimensional quantum dot in a perpendicular magnetic 
field B subject to a gate voltage V g . Unlike the traditional Coulomb blockade picture, it 
might contain several electron islands which can be regarded as metallic objects with induc- 
tive couplings among themselves. These are determined by a positive definite symmetric 
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matrix C whose diagonal elements Cu = C{ - the corresponding capacities - are positive, 
whereas the nondiagonal elements = Cjt, (i ^ j) - the corresponding coefficients of 
induction - are negative. The electrostatic energy of such a system may be written as 

E - = \ll Pv^i ~ V * N > (!) 

where iV, is the number of electrons on island z (the number of islands K > 1 might be 
around 10), iV = SiLi is the total numbers of electrons, and p = (pij) is the (symmetric 
positive-definite) matrix C . 

Beside the electrostatic energy it is assumed that electrons in each island occupy single 
particle quantum states (orbitals), whose energies €i a {i = 1,2, ..if; a = 1,2,..) depend on 
the confining potential as well as on the magnetic field. The latter is manifested through 
its orbital effects as well as due to Zeeman splitting (in which case the quantum number a 
contains also a spin label). The corresponding occupation numbers be either or 1. 

The system described above might then be represented by a classical Hamiltonian 

H = H c + H sp , (2) 

where the charging Hamiltonian H c is just the electrostatic energy of (|I|) written in terms 
of the orbital occupation numbers 

^ K K 

H c=oYl Pa E n i»] E n j*'} -^EE n «- ( 3 ) 

z i,j=l a a ' i=l a 

and the single particle part of the Hamiltonian H sp is 

K 

^sp ^ ] ^ ] ^ia^ia- (4) 
i=l a 

The precise form of the matrix elements = [C _1 ]jj as well as the single particle energies 
€i a are specified in the next section when we present our results. Despite the fact that the 
Hamiltonian H = H c + H sp is classical (and relatively simple), the elucidation of its spectrum 
for large iV and K is virtually hopeless. In order to compute the ground state energy E(N), 
one has to find the minimum of H on all the possible if -tuples (N\, N2, ■ ■ ■ , Nk) with the 
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constraint J2i=i = N. Note that the so-called "Coulomb Glass" model obtains as a 
special case when i refers to a lattice site with random energy e ia = e^, and a single orbital 
Ni = rii = 0, 1. The interaction matrix is then given by pij = 1/rij for i ^ j and pa = 0, 
where is the distance between sites i and j. 

The position of the N th conductance peak is given by the first difference of the ground state 
energy, namely, the chemical potential of the (isolated) dot, \i = E(N + 1) — E(N). The 
spacing between peaks is determined by the second difference (the inverse compressibility), 



The occurrence of close peaks for certain values of electron number N corresponds to small 
values of xn (recall that for a single island quantum dot, in which the single particle energies 
are neglected, the inverse compressibility is a constant 1/Cn). For the more general model 
described above the spacing distribution will of course fluctuate. In general, some constants 
appearing in the Hamiltonian H = H c + H sp are random (e.g. the elements of the matrix 
C and the single particle energies e^), but most experiments are performed on a single 
quantum dot, so that fluctuations are meant with respect to the electron number N. The 
numbers x(iV) might then be considered as values assumed by a random variable x which 
has a certain distribution function P(x)- 

It might be instructive to compute the distribution P(x) for a particular special case (albeit 
not realistic) which looks deceptively simple. The reader who is not interested may skip the 
rest of this section. Take K > 1, CV, = for i ^ j, and €i a = 0. The Hamiltonian H = H c is 
the classical energy of a set of K charged metallic bodies which are very far apart from each 
other, so there is no coupling among them. This problem then belongs to a class of problems 
dealing with statistics of spectra of independent quantum systems [[K]]. Here it refers to the 
addition spectra for a system composed of several independent subsystems, each of which 
has a certain spectrum which is a quadratic function of the number of particles residing on 
it. Thus, for example, by replacing C^ 1 — > hu>i, it can represent the energy of iV spinless 
fermions distributed among K independent oscillators. 




(5) 
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Our aim is to find the distribution function P(x)- ft charge is not quantized, the problem 
of finding the ground state energy is simply formulated as follows: Given K independent 
capacitors Cj, i = 1,...,K, (in general, with random capacitances Cj), how should the 
corresponding charges Qi be chosen so as to minimize the electrostatic energy E(Q) = 
Sj=i 2c" w hile holding the total charge Q = J2i=i Qi fixed? This constrained minimum 
problem is trivially solved by requiring dE(Q)/dQ = and expressing Q\= Q — Y J f=2 Qi- ft 
is evident from dimensional arguments that E(Q) is proportional to Q 2 , with some coefficient 
denoted by 1/2(7, in which C may be interpreted as the total capacitance of the system. 
Evidently, in this case the second derivative of E(Q) with respect to Q is independent of Q 
(and is easily calculated), namely d 2 E/dQ 2 = 1/(7. 

What happens then if charge is quantized? The formulation of the problem is now repeated 
albeit with Q = N,Qi = iVj (with N and iVj integers). Intuitively, one would expect P(x) to 
be centered about 1/(7, but, surprisingly, this is not the case, and P(x) is large near x — 0. 



An analytic expression for P(x) has been obtained recently by the authors [OJ]. In order to 
present it, let us assume that the K capacities (7j, % = 1, 2, . . . , K, are rearranged so that 
C\ > C 2 > C 3 ... > Ck- It is also useful to divide all the capacitors by the largest one, to 
obtain scaled capacitors q = GijG\. with 1 = c\ > c 2 > c 3 > . . . > ck- Now let us define 
the following quantities: 

a e = CiiCia—Cfc II C 1 _ c i)' ^ = 1,2, . . . ,iY - 1, (6) 

2<ii<i'2<..<ie<K j^im 



K 



c = E c - (7) 



1=1 



K 



B = H{l-a). (8) 

i=2 

With these notations we can provide a precise formula for the spacing distribution. The 
formula relates to the density function of the distribution and holds in the generic case 
(meaning that for a specific set of measure of parameters it is incorrect). The following is 



obtained if the spacings are normalized to lie in the interval [0, 1] (multiplying all of them 
by Cx): 

P (X) = ^f>^+ !)(! - XY" 1 + §5(1 ~ X). (9) 

A plot of P(x) 1S given in figure [I] and shows indeed that it has a maximum at x = and a 
delta function component at the inverse of the largest capacitance. In fact, when K —* oo, 
the weight of the delta function shrinks to zero and (after an appropriate renormalization) 
P(x) approaches the Poisson distribution e~ x . 



III. RESULTS 

We now return to the full Hamiltonian of Eq. |2|. It contains the coefficients of capacitance 
and induction matrix and the single particle energies as input. In choosing the actual 
numerical values we use a few guidelines, one of which is to avoid too many independent 
input data. As will be clear below, these are not fitting parameters, but rather a set of 
constants which are chosen once and for all on general physical grounds. First, for the 
electrostatic part, recall that the matrix C should be a symmetric positive definite matrix 
with Cu > and < for i ^ j. Besides, we expect it to be random. We then assume that 
Cjj are random numbers uniformly distributed between and W, whereas the nondiagonal 
elements are uniformly distributed between — w and 0. Actually, it is only the ratio w/W 
which matters, so one may assume W = 1 and use e 2 /W = 1 as an energy unit, leaving 
w as a single constant reflecting the strength of coupling between islands. Practically, for 
each realization of c^-'s drawn at random we check whether the resulting matrix is positive 
definite, and reject realizations which are not. Thus, we must put the restriction w < 1 
in order to generate a positive definite matrix C . Indeed, for w > 0.25 (and K = 5) most 
of the matrices generated randomly failed to be positive definite. Second, for the single 
particle energies, we consider each electron island i as a two-dimensional potential well 
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Vi(r) = ^Mufr 2 , where M is an effective mass. One may then regard cu i as a measure of 



the radius of the corresponding electron island. Since the capacitance Cu is also proportional 
to this radius we assume Ui = •yC^ 2 , where 7 is a constant reflecting the relation between 
charging energies and single particle energies. The single particle energies in each electron 
island (subject to a perpendicular magnetic field B) are then given explicitly. To be more 
specific, recall that in two dimensions there are two quantum numbers for the orbital motion 
(denoted hereafter by n, to) to which we add a spin index a — ±1. Then, with a = (n, to, a), 
we have 



where u c = eB/Mc is the cyclotron frequency, g is the g-factor and /xg is the Bohr magne- 
ton. Since g contains the effective mass it is not known accurately. Its value is constrained 
on physical grounds (see below). Note that in this scheme the spacings between single par- 
ticle energies are deterministic, and do not follow the Wigner surmise. The main cause of 
fluctuation is then due to the combination of non-random single particle energies and the 
occurrence of numerous charging energies. Finally, the number of electron islands K is de- 
termined by the size of the quantum dot. The four input data of the model are then K, w, 7 
and g. Note that the gate voltage V g does not have an important role here. Indeed, in 
actual experiments the variation of gate voltage serves to adjust the energies E(N) with the 
chemical potential of the leads, but here the ground state energies are calculated directly. 
In order to avoid redundancy we are content with having a single set of these constants 
which is physically reasonable. In particular, it assures that the charging energy is much 
larger than single particle level spacings and that the Zeeman splitting is small at moderate 
magnetic fields. Specifically, for B = we test the distributions for w = 0.03, 0.08 and 
0.20, while for B ^ we fix w = 0.03. Besides, we take 7 = 0.1, and g = 0.0081 is chosen 
such that for moderate fields the Zeeman splitting is of the order of the mean level spacing. 
Finally, the number of electron islands is fixed as K = 5. Indeed, we use just a brute force 
trial algorithm, and hence cannot treat systems with a large number of electron islands. 




(10) 
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With these prescriptions the ground state energies E(N) and their first and second differ- 
ences are calculated for B = up to N = 200 and for B ^ up to iV = 50. 
The first question we addressed is how the electrons are added among the islands. Figure |2| 
shows electron numbers Ni in each of the five islands, v.s the total electron number N, in the 
range 60 < N < 80 (for w = 0.03). Evidently, the order of curves is according to the value of 
the capacitance Cu. On a larger scale, the numbers Ni grow linearly with N, as it should be. 
Let us then consider the question of redistribution. From figure ^| we see that redistribution 
occurs only once, as N grows between 67 and 68 (see the vertical line). The addition of 
an electron to the second (or third) island involves also a transfer of an electron from the 
fourth island to the second (or the third) one. This scenario occurs also in other ranges of 
iV with the same proportion (namely about four percent). Thus, within the present model, 
redistribution is present, although it is rare and minimal. 

The next question we address is how the ratio between the inductive and capacitive coef- 
ficients (represented here by the parameter w) affects the distribution P(x)- It was shown 
in the previous section that for w = the distribution has a maximum at \ = 0. It is 
intuitively expected (based on the concept of level repulsion) that for w ^ the maximum 
will be shifted away from zero. This is also verified by our numerical results but, somewhat 
unexpectedly, the peak at small x persists even at w — 0.20. The distribution P(x) is drawn 
in figure |3] for w = 0.03,0.08 and w/W = 0.20. The normalization is such that the largest 
capacitor is Cu = 1.0, so that if other capacitors are absent the distribution would have a 
delta function at \ = 1. In all three cases the peak position is much smaller than 1. Recall, 
however, that besides the capacitance induction matrix C the total energy is determined also 
by the single particle energies in each island. The remarkable point is that the distribution 
is not Gaussian (which is the hallmark of spacing distributions in small quantum dots). 
Finally, we check the behavior of the first difference E(N + 1) — E(N) as a function of the 
magnetic field. Recall that this quantity is proportional to the position of the N th Coulomb 
blockade peak. As a measure of the strength of the magnetic field we use the parameter 
uj c /ujq, where uq is the harmonic oscillator frequency of the largest island. The positions 
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of the peaks for 39 < iV < 48 and w/W = 0.03 are displayed in figure |j. A comparison 
with the results displayed in figure 2 of Ref. J7| indicates a remarkable qualitative agree- 
ment. In particular, it shows that groups of two (and sometimes even three) electrons can 
tunnel through the quantum dot at almost the same gate voltage. The oscillations at small 
magnetic field just mark transitions to lower Landau levels as the magnetic field increases. 
The phenomena of alternate bunching (N,N + 1) — > (N — 1,N) is also reproduced in the 
present picture. 

In conclusion, we suggest a classical model in which a large semiconductor quantum dot 
is viewed as a collection of metallic electron islands with capacitive and inductive coupling 
among them. The effect of adding a magnetic field is manifested through its orbital as well 
as its spin contents. The model can explain the occasional occurrence of couples or even 
triples of closely spaced Coulomb blockade peaks, as well as the qualitative behavior of peak 
positions with the magnetic field. The results of the previous section displayed in figure [I] 
together with an analysis of the results displayed in figure |3| can provide an answer to the 
question why, in the present model, there are couples and sometimes even triples of close 
Coulomb blockade peaks. When the capacitors are completely independent, the distribu- 
tion is close to Poissonian as can be proved analytically. Introducing an inductive coupling 
spoils this picture, but does not destroy it completely, so that even at w = 0.2 there is still 
a probability for the occurrence of close peaks. The eventual decrease of the distribution 
near \ = seems to be related to the fact that, unlike the case of independent systems, 
the capacitance matrix is non-diagonal. Switching on the coupling here then has an effect 
similar to a weak "spacing repulsion" , similar to the familiar effect of perturbation in two 
level systems. 
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IV. FIGURE CAPTIONS 



FIG. 1. Distribution of inverse compressibility P(x) f° r a system containing independent subsys- 
tems each of which has a ground state energy proportional to the square of the number of particles 
it contains. The graph corresponds to a special case of equation |9] with K = 5. 

FIG. 2. Distribution of electrons among the five islands as a function of the total electron 
number N at zero magnetic field. The choice of constants K = 5, w = 0.03 is explained in the 
text. The number iVj of electrons on island i is commensurate with its capacitance Ca. In most 
cases the addition of an electron does not perturb the occupation of other islands. An example 
of redistribution is marked with a vertical line. An addition of an electron causes a minimal 
redistribution. 

FIG. 3. Effect of strength of inductive coupling w on the distribution P(x) of level spacings in 
the dot at zero magnetic field (not normalized). It is displayed for K = 5 and w = 0.03,0.08 and 
0.20 in (a), (b) and (c) respectively. 

FIG. 4. Coulomb blockade peak position for electron number N between 39 and 48 as functions 
of the magnetic field. Here uj is the oscillator frequency of the largest island and lo c is the 
cyclotron frequency. The units on the ordinate are energy units as explained previously. They are 
proportional to the gate voltage appropriate for the corresponding peak. The values chosen for the 
constants w = 0.03, K = 5, 7 = 0.1 and g = 0.0081 are explained in the text. 
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